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The halo model (HM) describes the inhomogeneous universe as a collection of halos. The full 
nonlinear power spectrum of the universe is well approximated by the HM, whose prediction can 
be easily computed without lengthy numerical simulations. This makes the HM a useful tool in 
cosmology. Here we explore the lensing properties of the HM by use of the stochastic gravitational 
lensing (sGL) method. We obtain for the case of point sources exact and simple integral expressions 
for the expected value and variance of the lensing convergence, which encode detailed information 
about the internal halo properties. In particular a wide array of observational biases can be easily 
incorporated and the dependence of lensing on cosmology is properly taken into account. This 
simple setup should be useful for a quick calculation of the power spectrum and the related lensing 
observables, which can play an important role in the extraction of cosmological parameters from, 
e.g., SNe observations. Finally, we discuss the probability distribution function of the HM which 
encodes more information than the first two moments and can more strongly constrain the large- 
scale structures of the universe. To check the accuracy of our modelling we compare our predictions 
to the results from the Millennium Simulation. 

PACS numbers: 98.62. Sb, 98.65.Dx, 98.80.Es 



I. INTRODUCTION 

Gravitational lensing by large-scale inhomogeneities 
affects the light from distant objects in a way which es- 
sentially depends on size and composition of the struc- 
tures through which the photons pass on their way from 
source to observer. The fundamental quantity describing 
this statistical (de)magnification is the lensing probabil- 
ity distribution function (PDF), which has to be under- 
stood well if one is to use cosmological observations to 
accurately map the expansion history and determine the 
precise composition of the universe (see, for example, [1- 
5]). It is not currently possible to extract the lensing 
PDF from the observational data and we have to resort 
to theoretical models. Two possible alternatives have 
been followed in the literature. A first approach (e.g. [6- 
9]) relates a "universal" form of the lensing PDF to the 
variance of the convergence, which in turn is fixed by the 
amplitude of the power spectrum, erg . Moreover the coef- 
ficients of the proposed PDF may be trained on some spe- 
cific N-body simulations. A second approach (e.g. [10- 
12]) is to build ab-initio a model for the inhomogeneous 
universe and directly compute the relative lensing PDF, 
usually through time-consuming ray-tracing techniques. 
The flexibility of this method is therefore penalized by 
the increased computational time. 

In Refs. [13, 14] we introduced the stochastic gravita- 
tional lensing (sGL) method which combines the flexi- 
bility in modeling with a fast performance in obtaining 
the lensing PDF. The sGL method is based on the weak 
lensing approximation and generating stochastic config- 
urations of inhomogeneities along the line of sight. A 
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numerical code based on sGL, as the publicly available 
turboGL package [15], can compute the lensing proba- 
bility distribution function for a given inhomogeneous 
model in a few seconds. The speed gain is actually a sine- 
qua-non for likelihood approaches, in which one needs 
to scan many thousands different models (see, e.g., [4]). 
This makes sGL a useful tool to study how lensing de- 
pends on cosmological parameters and how it impacts 
the observations. The method can also be used to simu- 
late the effect of a wide array of systematic biases on the 
observable PDF. 

Here we develop a self-consistent setup to easily calcu- 
late matter power spectrum and lensing properties of the 
desired model universe. The basic idea is to apply the 
sGL method to the so-called "halo model" (HM) (see, for 
example, [16-23]), where the inhomogeneous universe is 
approximated as a collection of different types of halos 
whose positions satisfy the linear power spectrum. We 
will give explicit integral expressions for mean and vari- 
ance of the lensing convergence, which will be directly 
related to the relative matter power spectrum. These 
results will be valid for point sources, while extended 
sources will be treated in forthcoming work. Our dis- 
cussion applies, for example, to the narrow light bundles 
emitted by distant Supernovae (SNe). We expect indeed 
that our formulas will be useful in the analysis of SNe 
observations as the convergence variance can be straight- 
forwardly included in the standard analyses. We also 
discuss the PDF of the HM which encodes more informa- 
tion than the mean and variance and can more strongly 
constrain the large-scale structures of the universe. 

This paper is organized as follows. In Section II we 
introduce the HM and the basic formalism to calculate 
its power spectrum, while in Section HI we will discuss 
its lensing properties. To check the accuracy of our mod- 
elling we compare in Section IV our predictions to the 
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results from the Millennium Simulation (MS) [24]. Fi- 
nally, we give our conclusions in Section V. 



II. THE HALO MODEL 

The halo model assumes that on small scales (large 
wavenumbers k) the statistics of matter correlations are 
dominated by the internal halo density profiles, while on 
large scales the halos are assumed to cluster according to 
linear theory. In other words, the nonlinear evolution is 
assumed to produce only concentrated halos. The model 
does not include intermediate low density structures such 
as filaments and walls. The two components are then 
combined together. We do it here by simple addition so 
that the total power spectrum is 



Pik,z)=PLik,z) + PH{k,z) 



(1) 



The first term on the right-hand side is also called 2- 
halo component and the second term 1-halo component. 
The usefulness of the halo model stems from the fact that 
both terms in Eq. (1) can be computed without having to 
resort to numerical simulations. We would like to stress 
that more sophisticated and accurate versions of the HM 
are available in the literature [25]. Here we use the sim- 
plest version of the HM as we merely wish to illustrate 
its good agreement with the full power spectrum. 



where the density Pmc = «o Pmo is the constant matter 
density in a co-moving volume. The halo function is by 
definition normalized to unity: J f{M, z)dM — 1 and 
we defined dn as the number density of halos in the mass 
range dM. The power spectrum due to randomly located 
halos is then [21] 



PH{k,z) 



dn{M, z) 



MWkiM, z) 



Pmc 



(4) 



where is the Fourier transform of the halo density 
profile: 



Wk{M, z) 



1 

M 
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p{r,M,z)— Anr^dr, (5) 



kr 



and R is the halo radius. For the halo profile p(r, M, z), 
we shall use the Navarro- Frenk- White (NFW) profile [29] . 
We will now apply this setup to the weak lensing. 



HI. LENSING 

The lens convergence k in the weak-lensing approxi- 
mation is given by the following integral evaluated along 
the unperturbed light path [30]: 



k{Zs) 



dr Pmc G{r,rs) 5M{r,t(r)) 



(6) 



Linear Power Spectrum 



The linear part of the power spectrum is as usual: 



where we defined the auxiliary function 
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Ho(^j T\klao)D\z), (2) 



where Hq and Oq are the present-day Hubble parameter 
and scale factor, rig is the spectral index and 5ho is the 
amplitude of perturbations on the horizon scale today, 
which we fix by setting a^. For the transfer function 
T{k) one can use the fits provided, e.g., by Ref. [2G]. Fit 
functions for the growth function D{z) could be found, 
e.g., in Refs. [27, 28], but it can also be easily obtained 
numerically [14]. 



Halo Power Spectrum 

The nonlinear or halo contribution has been obtained 
in, e.g., [18]. We will briefly summarize the important 
steps here. First we introduce the halo mass function 
f{M,z), which gives the fraction of the total mass in 
halos of mass M at the redshift z. The function f{M, z) 
is related to the number density n{M, z) by: 



dn(M, z) = n(M, z)dM 



Pmc 
M 



f{M,z)dM, (3) 



G(r,r,) 



fk{r)fk{rs - r) 
c^a fk{rs) 



(7) 



which gives the optical weight of an inhomogeneity 
at the comoving radius r. The functions a{t) and 
t{r) are the scale factor and geodesic time for the 
background FLRW model, and Ts = r{zs) is the co- 
moving position of the source at redshift z^. Also, 
/*;('') — sm{r\^)/y/k, r, smh{r^/—k)/^/—k depending 
on the curvature k >,=,< 0, respectively. Neglecting 
the second-order contribution of the shear, the shift in the 
distance modulus caused by lensing is expressed solely in 
terms of the convergence: 



Am{z) ~ 51ogio(l - k{z)) . 



(8) 



Eqs. (6) and (8) show that for a lower-than-FLRW col- 
umn density the light is demagnified (e.g., empty beam 
d]\i — —1), while in the opposite case it is magnified. 
In Eq. (6) the quantity SM{r,t) is the local matter den- 
sity contrast and it is clear that an accurate statistical 
modelling of the magnification PDF requires a detailed 
description of the inhomogeneous mass distribution. In 
this work we will model Sm according to the halo model. 
See [14] for a more refined modelling which also includes 
filamentary structures confining the halos. The following 
results will be valid for point sources (smoothing angle 
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6 = 0). We will discuss extended sources in a forthcom- 
ing paper. 

Because in Eq. (6) the contributions to the total con- 
vergence are combined additively, it is useful to decom- 
pose the density field into a sum of Fourier modes: 

SMir,t)= X (^e"'""'^M(k,i). (9) 

Eq. (9) can be used to separate the contributions to 
the convergence due to small-scale inhomogeneities (fc > 
kcut) from the ones due to large-scale inhomogeneities 
< fccut) so that, similarly to Eq. (1), it is: 

k{z) = KLiz) + kh{z) . (10) 

The quantity fccut separates the modes relevant to 
from the ones relevant to kh and can be defined as the 
scale at which the two power-spectrum components have 
the same value, Phikcnt, z) = Pnikcnt, z). 

We are in particular interested in the expected value 
and variance of the lensing convergence. From Eq. (10) 
it follows that: 

(k) = (kl) + (kh) , (11) 

<^1 = <^L+<^L- (12) 

We will now give directly computable expressions for 
Eqs. (11-12). 

Linear Contribution 

Because of photon conservation, if no lines of sight are 
obscured, the expected value of the weak lensing conver- 
gence is zero. This can be seen in Eq. (6) from the fact 
that {5M{r, t)) = 0. Mechanisms that are able to obscure 
lines of sight involve small scales and so we have 

(kl)=0. (13) 

The linear power spectrum, however, contributes to the 
variance of the convergence according to [31]: 

< = r drpl,aGHr,r,) ^P^(fc, .(.)) , (14) 

Jo Jo ^~ 

where the scale k^nt can depend on the redshift z. 

Halo Contribution 

The sGL method for computing the lensing conver- 
gence is based on generating stochastic configurations of 
halos and filaments along the line of sight and comput- 
ing the associated integral in Eq. (6) by binning into a 
number of independent lens planes. Because the halos are 
randomly placed their occupation numbers in parameter- 
space volume cells follow Poisson statistics. This allows 
rewriting Eq. (6) as a sum over these cells characterized 



by the corresponding Poisson occupation numbers. The 
various individual contributions to the convergence are 
then additively combined. By generating many halo con- 
figurations one can easily sample the convergence PDF. 
Besides the approximations relative to Eqs. (6-8), corre- 
lations in the halo positions are neglected here, similarly 
to Eq. (4) . In the full sGL one can model also filamentary 
structures confining the halos, thus accounting for some 
of the correlations among the halo positions. A detailed 
explanation of the sGL method can be found in [13, 14] 
and a publicly-available numerical implementation, the 
turboGL package, in [15]. 

In the present case the filamentary structures are 
turned off and the inhomogeneous matter distribution 
introduced in [14] exactly corresponds to the halo model. 
The sGL method then predicts the following direct and 
exact results for mean and variance of the convergence: 

{kh)= drG{r,r,) dn(M,0(r))x (15) 
Jo Jo 

l-R(M.z(r)) 

X / dA(6)(P,ur-l)S(^Af,z(r)), 

JO 

and 

<jI^= ArG\r,rs) dn(Af,z(r))x (16) 
Jo Jo 

X / dA(6)P,urS2(6,M,z(r)), 
Jo 

where the integral limits for the last two integrals are im- 
plicitly defined, dA{b) = 2Trbdb and S is the halo surface 
density: 

2rdr 

j:{b,M,z) = a^ , p(r,M,z). (17) 

Jb — 0^ 

Eqs. (15-16) can be directly integrated without having to 
consider the full formalism of the sGL method. We also 
would like to point out that, in contrast to, e.g., [32], our 
expressions use the halo profiles in real space and not 
in Fourier space, thus including higher order correlation 
terms beyond the power spectrum. As a direct conse- 
quence the survival probability Pgur can be included in 
the way the sGL method predicts. 

The quantity Pgui ~ Psui{b, M, z, Zg) is a generic func- 
tion that describes the probability that a light ray, which 
encounters a halo of mass M at the redshift z and im- 
pact parameter b for a source at redshift z^, does not 
fall below detection threshold. Pgui- can be used to simu- 
late the effect of a very wide array of systematic biases, 
such as any sources leading to obscuration of the light 
beam, either alone or in combination with restrictions 
arising from imperfect search efficiencies and strategies. 
For example, selection by extinction effects or by out- 
lier rejection mainly relates to high-magnification events 
which are clearly correlated with having high intervening 
mass concentrations (halos with large M and/or small b) 
along the light geodesic. Similarly, short duration events 
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might be missed by search telescopes or be rejected from 
the data due to poor quality of the light curve, for ex- 
ample in cases where a supernova is not separable from 
the image of a bright foreground galaxy. Probability of 
such events would be correlated with the brightness of 
the source, with the density and redshift of the interven- 
ing matter and of course with the search efficiencies. All 
these effects could be modelled by Pgur, carefully adjusted 
by the use of the astrophysical input. As a simple illus- 
tration, in the next Section we will model the survival 
probability by a step function in the impact parameter, 
such that the halo is opaque for radiuses smaller than the 
NFW scale radius Rs = R/c: 

where c(M, z) is the mass- and redshift-dependent NFW 
concentration parameter. 

Eqs. (15-16) allow to draw some general considerations. 
First, if Psur = 1, the expected value of the convergence 
is correctly zero as demanded by photon flux conserva- 
tion,^ showing the "benevolent" nature of weak lensing 
corrections for unbiased observations. If, however, the 
survival probability is not trivial, selection biases persist 
even in very large datasets so that the average conver- 
gence approaches a nonvanishing value. Second, Eq. (16) 
is a product of positive quantities and so a nontrivial 
survival probability always reduces the observed variance. 
Thus, neglecting existing systematic biases could lead to 
an underestimation of the observationally inferred vari- 
ance, and hence to too strong constraints on additional 
variance caused by inhomogeneities. 

Finally, Eqs. (15-16) are integrated (see second dn- 
integral) over all halo masses < M < oo. It is, how- 
ever, straightforward to generalize them by changing the 
integration limits in order to have expected value and 
variance relative to a specific halo mass bin AM. This 
may be useful if one wants to connect the lensing given by 
halos within a mass bin (see Fig. 1) to the corresponding 
amount of correlation. 



IV. COMPARISON WITH THE MILLENNIUM 
SIMULATION 

We will now compare the halo model power spectrum 
and lensing properties to those from the Millennium Sim- 
ulation. Accordingly, we will fix the cosmological param- 
eters to h = 0.73, = 0.25, Qb = 0.045, fl/^ = 0.75, 
w = —1, as = 0.9 and rig = 1; see [24] for more details. 
Let us point out that in [14] we have already accurately 
reproduced the MS lensing PDF with an sGL modelling 



This should be approximately true [3.3] if strong lensing events 
leading to secondary images and caustics [31] do not play a signif- 
icant role as far as the full-sky average is concerned [10, 35, 36]. 
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FIG. 1. Shown is tlie variance per mass bin for a source 
at redshift z = 1. The variance is computed using Eq. (16) 
restricted to the mass bin AAf =]10", 10"+^] H-'^Mq. In the 
plot n labels the halo mass according to M„ = lQ"h~^ Mq. 
For halos of mass M < IO^/i'^Mq the contribution to the 
total variance is negligible. See Section IV for more details. 



which included not only the halos, but also the filamen- 
tary structures. Our goal here is different; we are deliber- 
ately using a simpler modelling of the inhomogeneities in 
order to study the validity of the halo model for lensing 
analyses. 

Because we are assuming that all (virialized) matter is 
concentrated in halos, we need to use a halo mass func- 
tion whose integral is normalizable to unity. We therefore 
adopt the halo function provided by Sheth & Tormcn in 
Ref. [37], which should be approximately valid for the 
full mass range of the integrals in Eqs. (15-16). Halo 
functions obtained through numerical simulations, albeit 
possibly more precise, are valid only above a mass cut- 
off imposed by the numerical resolution of the simulation 
itself [38]. We point out, however, that halos with mass 
smaller than M^ut ^ IQ'^hr^ Mq act effectively as a mean 
field in weak lensing, and so mass functions from numeri- 
cal simulations which are valid down to Mcut may be used 
within the present setup by introducing explicitly Mcut 
in Eqs. (15-16). This fact can be appreciated quantita- 
tively by computing the convergence variance (16) per 
mass bin AM as shown in Fig. 1. Clearly for M < Mcut 
the contribution to the total variance is negligible which 
means that the corresponding halos behave effectively as 
a fine homogeneous dust as far as weak lensing is con- 
cerned [14]. 

In Fig. 2 we have plotted the power spectrum of 
Eq. (1), together with the halo and linear components. 
Also plotted is the power spectrum of the Millennium 
Simulation. The agreement is rather good, although the 
halo model power spectrum is slightly underpowered on 
small scales. A better power spectrum could be obtained 
with a more sophisticated HM, but this would not be rel- 
evant for the main goal of this paper, which is to propose 
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FIG. 2. Shown is the present-day power spectrum of the halo 
model (orange solid line and Eq. (1)) for the ACDM universe 
of the Millennium Simulation, together with the halo (blue 
dashed line and Eq. (4)) and linear (green dot-dashed line and 
Eq. (2)) components. Note how the halo component gives a 
constant contribution at large scales (shot noise). Also shown 
for comparison is the power spectrum (black dotted line) rel- 
ative to the Millennium Simulation [24]. 
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FIG. 3. Shown is the redshift dependence of the dispersions 
due to the linear (green dot-dashed line and Eq. (14)) and halo 
(blue dashed line and Eq. (16)) components, for the ACDM 
model of the Millennium Simulation. Also plotted is the sum 
(orange solid line) of the variances and the dispersion (black 
dotted line) relative to the Millennium Simulation [35, 36]. 
The labeling is as in Fig. 2. 



using Eqs. (15-16) for lensing. The oscillations at large 
scales are absent in the HM spectrum because we used a 
transfer-function fit that reproduces the baryon-induced 
suppression on the intermediate scales but ignores the 
acoustic oscillations, which are again not relevant for us 
here [26]. For the NFW concentration parameters we 
used the fit^ of Ref. [39]. Note also that, as said before, 



^ We increased the fit by a factor of 4.67/3.93 (see [39]) in order to 
match the results of [40] which are relative to the MS cosmology. 



it is straightforward to calculate the correlation per mass 
bin by binning Eq. (4) in the same way Eq. (16) has been 
binned in Fig. 1. 

In Fig. 3 we have plotted the redshift dependence of 
the convergence dispersions due to the linear and halo 
components, computed with Eq. (14) and Eq. (16), and 
also the sum of the variances which shows that the lin- 
ear contribution is practically negligible. We wish to 
stress that here we are studying the case of point sources, 
i.e., we are adopting a smoothing angle = 0. For ex- 
tended sources and/or finite smoothing angles the re- 
sults of Fig. 3 change, and for smoothing angles larger 
than some value 6 the linear contribution dominates over 
the one due to the halos. At present we cannot treat 
nonzero smoothing angles within the sGL model, but 
we can nonetheless estimate the order of magnitude of 
9. Assuming for simplicity a flat universe, the optical 
weight function G of Eq. (7) peaks at half the comov- 
ing distance to the source and vanishes at the observer 
and source locations. The physics at rs/2 is therefore 
the most relevant for lensing, and a smoothing angle 9 
can be converted into a smoothing scale by means of 
A = 9 dA{rs/2), where cIa is the angular diameter dis- 
tance. Let us now take, for example, a source at z = I 
and a smoothing scale of A ~ 1 Mpc, which should be 
enough to smooth out the inhomogeneous halo profiles. 
For the MS cosmology the corresponding smoothing an- 
gle is ^ ~ 3'. Also plotted in Fig. 3 is the dispersion 
relative'^ to the Millennium Simulation [35, 36], which 
shows a good agreement with the theoretical result. We 
wish to stress that it is very easy to evaluate the integrals 
of Eqs. (15-16) numerically, and their predictions can be 
straightforwardly implemented in analyses based on 
Gaussian likelihoods. However, it should be remembered 
that Gaussian analysis is justified if the lensing PDF is 
not strongly skewed. When this is not the case the full 
PDF has to be used for the likelihood analysis, as was 
done, e.g., in [4]. 

In Fig. 4 we have plotted the lensing PDFs obtained 
with the sGL method (histograms). Only the (l-)halo 
contribution is included as this is the direct output of 
the sGL method. Indeed at present we cannot properly 
include the linear contribution which, in any case, should 
be subdominant as shown in Fig. 3. For z = 0.83 and 
1.08 the agreement with the MS results (dotted lines) 
from [35, 36] is remarkable, especially considered the sim- 
plicity of the setup. For z = 1.50 the halo model predicts 
a less skewed PDF. This could mean that the model we 
are using to describe the halos (Sheth & Tormen mass 



In order not to exceed the validity region of the weak-lensing 
approximation, the dispersion relative to the MS has been calcu- 
lated in the weak-lensing regime k <S 1 by applying a cut at large 
magnifications. More precisely, to ensure that the PDF is basi- 
cally unaltered, the cut has been determined by asking that the 
high-magnification tail has the total probability of 10"'^. This 
resulted in Kcut = 0.063,0.127,0.195,0.255,0.365 for the PDFs 
relative to ^ = 0.28,0.56,0.83, 1.08, 1.50, respectively. 
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FIG. 4. Shown are the lensing PDFs for a source at (from 
taller to shorter graph) z = 0.83, 1.08 and 1.50 for the ACDM 
model of the Millennium Simulation. The dotted lines are 
the lensing PDFs generated by shooting rays through the MS 
[35, 36], while the histograms are the corresponding lensing 
PDFs obtained with the sGL method [1-3-15]. 
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FIG. 5. Shown is the redshift dependence of the conver- 
gence dispersion with (red dot-dashed line) and without (blue 
dashed line) the selection effects of Eq. (18), for the ACDM 
model of the Millennium Simulation. Also plotted is the mean 
convergence (red solid line) with selection effects (without se- 
lection effects it is zero). 



function and NFW concentration parameters from [39]) 
is less accurate at high redshifts. For example, a better 
agreement with [35, 36] is found if the concentration pa- 
rameters are increased by 40%, and similar results may 
be found by slightly changing the shape of the mass func- 
tion. Alternatively the extra skewness could be caused 
by unvirialized low-density large-scale structures like fil- 
aments and walls. These could be particularly important 
at high redshifts where lesser amounts of virialized struc- 
tures, which are the only ones accounted for by the halo 
model, are present. Indeed, a very good agreement with 
MS was found in [14] where the halo model was extended 
to include filamentary structures. Note, finally, that the 
differences in shape between the two PDFs at 2 = 1.50 
do not seem to strongly impact the dispersions plotted 
in Fig. 3. 



Selection effects 
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FIG. 6. Shown is the lensing PDF for a source at z = 1.08 
with (taller histogram) and without (shorter histogram) the 
selection effects of Eq. (18), for the ACDM model of the Mil- 
lennium Simulation. 



It is interesting to discuss the impact on lensing of the 
toy survival probability of Eq. (18). In Fig. 5 we show 
the redshift dependence of the convergence dispersion 
(16) with (dot-dashed line) and without (dashed line) 
the selection effect, which halves the resulting disper- 
sion. Also plotted is the (negative) mean convergence 
(15), which shows the possible bias present even in very 
large datasets. 

In Fig. 6 we show the full lensing PDF with selection 
effects included (taller PDF). The lensing PDF with uni- 
tary survival probability is also shown for comparison 
(shorter PDF). As expected the selection effects reduce 
the high magnification tail of the PDF, without sizeably 
changing the tail at low magnifications. See [14] for an- 



other discussion of the impact of selection effects on lens- 
ing. 

Finally, we would like to point out that a finite smooth- 
ing angle reduces the skewness of the PDF as it reduces 
the contrast in the matter surface density. If observa- 
tional biases are added, however, the overall effect on the 
PDF is more involved and depends strongly on the kind 
of selection effect considered. In the case of obscuration, 
for example, the bias on the mean convergence should 
remain basically unchanged as photons relative to un- 
seen lines of sight are lost no matter the smoothing angle 
adopted. Therefore for increasing the PDF should tend 
to a gaussian with decreasing dispersion peaked at the 
mean biased convergence. 
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V. CONCLUSIONS 

By applying the sGL method to the halo model, we 
have obtained a simple setup that allows to quickly com- 
pute power spectrum and lensing observables for a de- 
sired model universe. We wish to stress that there are no 
extra free parameters in the model, besides the ones rel- 
ative to the halo model derived from the literature (i.e., 
halo mass and profile fitting functions). In particular, 
we have given exact and directly-computable expressions 
for the expected value (see Eq. (15)) and variance (see 
Eq. (16)) of the convergence. We stress that Eqs. (15- 
16) contain more information than analogous expressions 
based on the power spectrum, like Eq. (14). They are a 
direct consequence of the sGL method which models the 
halo profiles in real space, and therefore include infor- 
mation of coherent structures described by higher order 
correlation terms beyond the power spectrum. In par- 
ticular, a wide array of systematic biases, which could 
be relevant for the extraction of cosmological parameters 
from, e.g., SNe observations and which persist even in 
very large data sets, can be included in Eqs. (15-16). We 
have shown a quantitative example of the impact of se- 
lection effects on lensing with Figs. 5 and 6. Moreover, 
these equations have the crucial advantage of including 
the dependence of lensing on cosmology, which was shown 
in [4] to give sizeable corrections to the confidence level 
contours from SNe observations. This is to be contrasted 
with the usual approach (e.g., [5, 41]) where lensing ef- 
fects are included by means of a cosmology-independent 
dispersion, which neglects the important dependence on, 
e.g., as, ^Im,a and h. 

We have then checked the accuracy of this setup 
against the cosmology of the Millennium Simulation. We 
have found that the theoretical predictions of the vari- 
ance are in good agreement with the results relative to 
the Millennium Simulation [35, 36]. This is most impor- 
tant as we wish here to propose to the community the use 
of Eqs. (15-16) for lensing, which can be easily included 
in the standard analyses. We have also found that the 
contribution of the linear (2-halo) component to the to- 



tal variance is practically negligible. Then we have shown 
the lensing PDFs of the halo model as computed using 
the sGL method and the turboGL package [15] for the 
redshifts of z = 0.83, 1.08 and 1.50. Given the simplic- 
ity of the setup, we found a remarkable agreement with 
the results from the Millennium Simulation, especially 
dX z = 0.83 and 1.08. One could conclude that, as far 
as the weak lensing magnification is concerned, the main 
contribution comes from the smooth halo profiles, with a 
smaller correction due to correlations among the halo po- 
sitions and substructures within the halos (even though 
one should keep in mind the possibility that the MS itself 
might lack some interesting substructures due to its res- 
olution limit). At z = 1.50 we have found that the halo 
model predicts a less-skewed PDF. This could be either 
due to the fact that the model we are using to describe 
the halos (Sheth & Tormen mass function with NFW ha- 
los) is less accurate at high redshifts or due to unvirial- 
ized low-density large-scale structures like filaments and 
walls, which could be important at higher redshifts where 
less virialized structures are present. One could indeed 
improve the modelling of the inhomogeneities, for exam- 
ple by explicitly introducing the filamentary structures as 
in [14], or by considering substructures within halos and 
filaments. Indeed, any extra matter structures such as fil- 
aments and walls increase the skewness of the PDF, and 
it would be interesting to see to what extent one can al- 
ter the lensing PDF by introducing such structures with- 
out conflicting with the observational constraints on the 
matter power spectrum. We will develop these thoughts 
further in a forthcoming paper. 
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